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In  what  follows,  the  results  of  the  research  under  this  program  are  briefly  summarized. 
More  details  can  be  found  in  the  preprints  of  submitted  papers  that  are  attached  at  the  end  of  this 
report. 

1.  Generation  of  extended  and  dense  plasma  channels  in  air  using  igniter-heater 
excitation  scheme 

Femtosecond  laser  filamentation  allows  for  a  substantial  degree  of  control  over  plasma 
generation  in  the  ambient  air.  The  longitudinal  position  of  the  plasma  filament  can  be 
straightforwardly  manipulated  by  introducing  a  negative  temporal  chirp  onto  the  laser  pulse.  As 
the  chirped  pulse  propagates  in  the  dispersive  air,  it  is  gradually  compressed  in  time  domain, 
reaching  its  maximum  peak  power  at  a  particular  propagation  distance  determined  by  the  value  of 
the  chirp.  The  placement  of  the  plasma  channels  within  the  transverse  beam  profile  can  be 
controlled  by  beam  shaping. 

Longitudinally  extended  plasma  channels  generated  in  air  through  femtosecond  laser 
filamentation  have  numerous  potential  applications  such  as  the  guidance  of  electricity  and 
emission  of  THz  radiation.  However,  the  density  of  free  electrons  generated  through  femtosecond 
filamentation  in  air,  under  typical  conditions,  is  only  about  1022  m'3,  which  is  too  low  for  many 
applications. 

In  this  program,  we  have  conducted  preliminary  experiments  on  the  generation  of  dense 
plasma  channels  through  the  combined  application  of  ultraintense  femtosecond  and  high-energy 
nanosecond  laser  pulses.  The  laser  hardware  used  in  this  research  has  been  acquired  under  the 
support  of  two  rounds  of  AFOSR  DURIP  program,  and  the  theoretical  support  has  been  provided 
through  the  interaction  with  the  AFOSR-supported  MURI  on  laser  filamentation. 

The  experimental  setup  used  in  this  work  is  schematically  shown  in  Figure  1.  The  optical 
breakdown  of  the  ambient  air  is  driven  by  10  ns-long  heater  pulses  generated  by  a  Q-switched 
Nd:YAG  laser  operating  at  1064  nm  wavelength.  The  maximum  pulse  energy  of  the  heater  laser 
is  3.3  J.  The  heater  pulses  are  combined  with  femtosecond  igniter  pulses  at  800  nm  wavelength. 
The  igniter  pulses  are  temporally  chirped  to  about  500  fs  duration  in  order  to  avoid  filamentation 
and  associated  damage  in  the  dichroic  beam  combiner.  The  heater  beam  is  focused  by  an  axicon 
lens  with  the  apex  angle  of  175°.  The  femtosecond  igniter  is  weakly  focused  by  an  ordinary  lens 
with  the  focal  length  of  0.5  m.  Prior  to  the  lens,  the  igniter  beam  is  apertured  in  order  to  clear  a 
hole  made  on  the  axis  of  the  axicon  lens  used  to  focus  the  heater. 


Femtosecond  igniter  pulse, 
600 nm  wavelength, 
lOmJ  per  pulse 


Figure  1:  Schematic  of  the  igniter-heater  setup  for  dense  plasma 
channel  generation  in  air. 


Using  this  setup,  we  have  demonstrated  that  the  density  of  the  dilute  plasma  left  in  the  wake 
of  the  femtosecond  laser  filament  can  be  boosted  by  a  factor  of  up  to  200  by  the  application  of  the 
properly  timed  nanosecond  pulse  with  Joule-level  pulse  energy.  Quite  naturally,  the  amount  of  the 
plasma-density  enhancement  depends  on  the  delay  between  the  two  pulses.  The  enhancement 
shows  clear  threshold-like  dependencies  on  the  energies  of  both  heater  and  igniter  pulses. 

This  behavior  is  currently  being  modeled  by  the  participants  of  the  AFOSR-supported  MURI 
on  laser  filamentation.  The  preliminary  modeling  results  indicate  that  heating  of  the  electron 
plasma  by  the  heater  pulse  plays  a  major  role  in  the  development  of  the  optical  gas  breakdown, 
as  the  electron  temperature  is  directly  related  to  the  recombination  rate  of  the  plasma.  We 
anticipate  that  the  effect  will  be  dramatically  different  in  an  atomic  gas  such  as  Argon.  We  have 
conducted  preliminary  experiments  on  the  igniter-heater  technique  in  Argon  atmosphere  and 
indeed  observed  a  different  breakdown  regime  from  that  in  the  air.  This  work  is  currently  in 
progress.  Our  ultimate  goal  is  to  develop  the  igniter-heater  approach  into  a  viable  technique  for 
the  controllable  generation  of  dense  plasma  channels  and  channel  arrays  that  would  be  suitable 
for  such  applications  as  guiding  microwaves. 


2.  Filamentation  of  femtosecond  Bessel  beams  of  higher  order 

Early  on  in  filamentation  studies,  it  was  realized  that  beam  shaping  offers  an  effective  means 
of  control  over  the  extent  of  the  generated  filaments  and  their  placement  within  the  beam  profile. 
The  application  of  various  beam  shapes  to  filament  formation,  in  both  gaseous  and  condensed 
media,  has  been  extensively  studied  in  the  past.  Examples  of  beams  previously  used  for  the 
formation  of  filaments  with  particular  properties  include  fundamental  Bessel  beams,  vortex 
beams,  necklace  beams,  and  self-bending  Airy  beams. 

We  have  conducted  experiments  on  filamentation  of  femtosecond  Bessel  beams  of  higher 
order  in  water.  In  addition  to  being  optical  vortices,  these  beams  are  diffraction-resistant.  They 
have  bottle-like  transverse  intensity  distributions  that  retain  their  shapes  on  propagation  over 
extended  distances.  The  goal  of  our  study  is  to  use  these  properties  of  higher-order  Bessel 
beams  for  the  generation  of  extended,  bottle-like  distributions  of  filaments  that  may  be  useful  in 
applications. 


J4 

CO 

“ > 

Ji2 

c_ 

OO 

0.26  mJ 

• 

0.32  mJ 

A 

0.37  mJ 

• 

•  • 

•  • 

0.6  mJ 

.  •  *  • 

• 

• 

% 

%  0 

200jtm 

%  # 

•  # 

'  t 

•  » 

•  • 

f 

• 

*  •  «•  * 

0.76  mJ 

0.77  mJ 

0.88  mJ 

1.1  mJ 

i . v 

v 

••  7. 

o 

v 

Figure  2:  Filamentation  of  Bessel  beams  of  various  orders  in  water.  The 
beams  are  generated  from  a  fundamental  Gaussian  beam  by  a 
programmable  phase  mask. 

Examples  of  the  filament  patterns  formed  by  Bessel  beams  of  various  orders  are  shown  in 
Figure  2.  The  self-focusing  collapse  is  accompanied  by  the  emission  of  the  forward-propagating 
white-light  continuum  radiation. 


If  the  dominant  rings  of  the  Bessel  beams  were  perfectly  smooth,  2n+1  filaments  would  be 
expected  to  form  on  the  main  ring  of  the  nth-order  beam  due  to  the  azimuthal  modulation 
instability.  However,  in  our  case,  the  placement  of  filaments  is  forced  by  the  beaded  structure  of 
the  beams,  which  is  caused  by  the  27r-modulo  phase  mask  that  we  use  to  generate  these  beams. 
As  the  pulse  energy  is  increased  further,  filaments  start  forming  on  the  secondary  rings  of  the 
beams,  as  shown  in  the  bottom  row  of  Figure  2. 

Our  experiments  in  water  involved  femtosecond  Bessel  beams  with  moderate  energy  ~1  mJ. 
We  are  presently  extending  this  work  onto  filamentation  of  much  more  energetic  Bessel  beams  in 
air.  To  create  higher-order  Bessel  beams  with  high  energy,  we  use  static  phase  masks  fabricated 
for  us  by  an  affiliated  holography  group  at  the  University  of  Arizona.  The  involvement  of  that 
group  into  this  work  was  partially  supported  through  this  program  which  paid  for  materials  and 
supplies  as  well  as  for  the  part-time  support  of  the  technician  involved  in  the  mask  fabrication. 
The  major  challenge  of  the  work  on  filamentation  of  higher-order  Bessel  beams  in  air  is 
associated  with  the  resistance  of  these  beams  to  filamentation.  The  total  energy  of  femtosecond 
laser  pulses  available  to  us  is  limited  to  about  35  mJ,  and  only  few  percent  of  that  energy  is 
carried  by  the  center  ring  of  a  higher-order  Bessel  beam  where  the  filaments  are  formed.  This 
work  is  presently  on-going. 


3.  Generation  of  sophisticated  spatio-temporal  pulse  waveforms  including  various 
light  bullets  and  application  of  these  waveforms  to  filamentation  studies. 

In  order  to  develop  techniques  for  control  over  filament  formation  and  dynamics,  we  have 
recently  started  experimenting  with  temporal  shaping  of  femtosecond  laser  pulses.  The  particular 
pulse  shape  that  we  have  studied  in  detail  is  the  accelerating  Airy  pulse.  Our  initial  experiments  in 
this  area  focused  on  the  generation  of  these  pulses,  their  characterization,  and  application  of 
these  pulses  to  supercontinuum  generation  in  optical  fibers.  This  is  a  relatively  simple  one¬ 
dimensional  problem,  in  which  the  transverse  beam  profile  is  invariant  on  propagation  as  it  is 
defined  by  the  properties  of  the  optical  fiber  where  lasers  pulses  propagate. 
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Figure  3:  Top  row,  left:  Temporal  waveforms  of  the  Airy  pulse  with 
-60,000  fs3  cubic  spectral  phase  after  6  cm  propagation  in  the  fiber.  The 
main  peak  of  the  Airy  waveform  has  produced  a  soliton-dispersive  wave 
pair  and  is  nearly  regenerated.  Right:  Same  after  16  cm  propagation.  The 
regenerated  main  peak  of  the  pulse  has  produced  another  soliton- 
dispersive  wave  pair.  Bottom  row:  Corresponding  spectra. 


We  have  conducted  both  numerical  simulations  and  experiments  on  the  propagation  of  Airy 
pulses  in  a  highly-nonlinear  optical  fiber  and  found  that  the  pulse  propagation  dynamic  is  different 


for  the  cases  when  the  tail  of  the  Airy  waveform  propagates  in  front  of  and  behind  the  main  peak 
of  the  pulse.  In  the  first  case,  the  dominant  peak  of  the  pulse  recreates  itself  after  several  soliton 
formation  events,  resulting  in  the  generation  of  distinct  dispersive-wave  features  in  the  output 
spectrum.  In  the  second  case,  the  soliton  generated  from  the  first  collapse  of  the  main  peak  of  the 
pulse  passes  through  the  tail  of  the  waveform  and  disrupts  its  ability  to  self-regenerate.  The 
results  of  the  modeling  are  in  qualitative  agreement  with  the  experiments.  The  temporal  dynamics 
and  supercontinuum  spectra  generated  by  an  Airy  pulse  propagating  tail-first  are  shown  in  Figure 
3. 

We  have  recently  extended  this  work  on  the  studies  of  filamentation  with  self-healing  Airy 
waveforms  in  water.  The  examples  of  the  angularly-resolved  spectra,  for  the  cases  of  Airy  pulses 
propagating  tail-first  and  peak-first  are  shown  in  Figure  4.  These  patterns,  in  which  the  emission 
wavelength  is  resolved  along  the  horizontal  axis  and  the  emission  angle  with  respect  to  the 
geometrical  beam  axis  is  resolved  along  the  vertical  axis,  are  indicative  of  the  complex  temporal 
evolutions  of  these  pulse  waveforms  on  propagation  in  the  water.  Similar  to  the  case  of  the 
supercontinuum  generation  with  Airy  pulses,  the  orientation  of  the  tail  of  the  pulse  affects  the 
generated  spectral  features.  We  are  presently  conducting  extensive  numerical  simulations  of  this 
complex  propagation  regime. 
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Figure  4:  Angularly  resolved  spectra  of  the  broadband  forward- 
propagating  emission  from  filaments  generated  by  femtosecond  temporal 
Airy  pulses  in  water.  The  pulse  energies  and  the  values  of  the  cubic 
phase  are  indicated  in  the  individual  figures.  The  case  of  a  Gaussian 
pulse  is  shown  in  the  left  panel.  This  pattern  is  indicative  of  a  single  pulse¬ 
splitting  event  that  accompanies  filamentation  in  this  case.  In  the  case  of 
an  Airy  pulse  with  the  tail  propagating  first  (center  panel),  there  are  two 
pulse-splitting  events  that  result  from  the  self-healing  nature  of  the  pulse 
waveform.  In  the  case  of  the  same  waveform  but  with  the  tail  propagating 
behind  the  main  peak,  the  second  splitting  event  is  smeared  out  due  to 
the  interference  between  the  splitting  products  and  the  tail  of  the 
waveform.  This  interference  inhibits  the  pulse  self-regenerating  ability. 


4.  Generation  of  harmonics  of  tightly  focused  mid-IR  femtosecond  pulses  in  ambient 
air. 

This  work  targeted  the  implementation  of  the  proposal  put  forth  by  the  AFOSR-supported 
MURI  team  back  in  2010.  The  goal  was  to  validate  or  disprove  the  recently  suggested  high-order 
Kerr  effect  theory  that  claims  that  the  conventional  Kerr  effect  changes  sign  at  high  light 
intensities  of  the  order  of  30  TW/cm2.  If  that  suggestion  is  proven  valid,  then  the  generation  of  3rd 
and  5th  harmonic  radiation  that  accompanies  laser  filamentation  would  become  very  efficient.  Our 
results  indicate  that  this  is  not  the  case,  in  agreement  with  the  established  paradigm  of  ultrafast 
nonlinear  optics  that  is  based  on  the  conventional  Kerr  effect. 


grating 

Figure  5:  Experimental  setup  for  3rd  and  5th  harmonic  generation  in 
ambient  air  using  mid-IR  pulses  at  2.2pm  wavelength,  in  a  tightly  focused 
geometry. 

Our  experimental  setup  is  schematically  shown  in  Figure  5.  Femtosecond  pulses  at  2.2  pm 
center  wavelength  are  produced  by  a  commercial  high-energy  Optical  Parametric  Amplifier  (OPA, 
Palitra-HE  model  by  Quantronix  Corporation).  The  OPA  is  pumped  by  about  40  fs-long,  20  mJ 
pulses  at  800  nm  wavelength  generated  by  an  ultrafast  Ti:Sapphire  amplifier  chain  operating  at  a 
10  Hz  pulse  repetition  frequency.  The  OPA  is  tuned  to  produce  signal  and  idler  pulses  at  1.26  pm 
and  2.20  pm  center  wavelengths,  respectively.  The  output  beam  from  the  OPA  has  a  diameter  of 
about  10  mm.  Both  signal  and  idler  pulses  have  energies  of  about  1.5  mJ.  The  duration  of  the 
pulses  is  estimated  by  the  OPA  manufacturer  at  between  40  and  60  fs. 

The  third  and  fifth  harmonics  of  the  2.2  pm  pump  are  generated  in  the  ambient  air, 
near  the  focal  plane  of  the  focusing  lens.  The  estimated  intensity  of  the  pump  beam  near  the 
focus  is  above  the  turnover  point  for  the  alleged  sign  reversal  of  the  Kerr  effect  in  the  HOKE 
model  and  also  above  the  ionization  threshold  of  air.  The  length  of  the  short  and  faint  plasma 
spark  generated  in  the  vicinity  of  the  interaction  zone  is  less  than  1  cm.  Under  these  conditions, 
the  ratio  of  the  third  and  fifth  harmonic  yields,  according  to  the  high  order  Kerr  effect  theory, 
should  approach  one. 

After  collimation  with  another  lens,  the  signals  at  third  and  fifth  harmonics  are  spatially 
separated  by  a  300  lines  per  millimeter  transmissive  diffraction  grating.  The  spectra  of  the  two 
harmonics  are  recorded  by  a  fiber-coupled  spectrometer  (model  USB4000  by  Ocean  Optics 
Corporation),  with  the  aid  of  two  additional  lenses  that  separately  focus  the  harmonic  beams  after 
they  are  separated  by  the  grating.  These  spectra  are  shown  in  Figure  6.  They  are  centered  at  733 
nm  and  438  nm,  very  close  to  one-third  and  one-fifth  of  the  pump  wavelength  of  2.20  pm. 
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Figure  6:  Measured  spectra  of  the  third  and  fifth  harmonics. 

The  energies  of  the  two  harmonic  pulses  are  measured  by  dedicated  silicon  detectors  that 
were  calibrated  against  a  pyroelectric  detector,  using  reference  laser  pulses  at  800  nm 
wavelength.  Accounting  for  all  losses  experienced  by  the  generated  harmonics  in  our  setup,  the 
energies  of  the  third  and  fifth  harmonics  are  found  equal  to  1.9  pJ  and  0.4  nJ,  respectively.  These 
pulse  energies  correspond  to  yields  of  about  1.5  x  10'3  for  third  harmonic  and  3.0  x  10'*  for  fifth 
harmonic,  with  respect  to  the  energy  of  the  pump  pulse  at  2.2  pm  wavelength  in  the  interaction 
zone.  The  ratio  of  yields  for  the  two  harmonics  is  about  2x10^,  in  contradiction  with  the  higher- 
order  Kerr  effect  theory. 
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Experiments  on  gas  breakdown  by  powerful  optical 
pulses  date  back  to  the  early  days  of  lasers.  Shortly 
after  the  demonstration  of  Q-switching  of  a  solid-state 
laser  in  the  early  1960’s  [1],  a  Q-switched  ruby  laser  was 
used  for  the  generation  of  plasma  sparks  in  various  noble 
gases  [2] .  Subsequent  studies  of  laser  breakdown  in  gases 
were  primarily  motivated  by  the  potential  use  of  this  phe¬ 
nomenon  in  controlled  nuclear  fusion.  (For  a  review  of 
the  subject  see,  for  example,  [3].) 

The  mechanisms  responsible  for  gas  breakdown  by 
nanosecond  laser  pulses  are  well  understood.  During 
breakdown,  free  electrons  are  accelerated  by  the  laser 
field  through  the  inverse  bremsstrahlung  process  which 
involves  their  collision  with  neutral  and  already  ionized 
molecules.  Once  the  accelerating  electrons  attain  en¬ 
ergy  that  equals  the  ionization  potential  of  the  neutral 
molecules,  these  electrons  impact-ionize  the  molecules 
forming  an  electron  avalanche.  The  few  seed  electrons 
needed  to  start  the  avalanche  are  typically  provided  by 
the  multi-photon  ionization  of  the  easily  ionizable  impu¬ 
rities  in  the  gas. 

At  a  given  pulse  duration,  optical  breakdown  by  a 
nanosecond  laser  pulse  is  a  threshold-like  process  with 
respect  to  the  peak  intensity  of  the  pulse.  The  numeri¬ 
cal  value  of  this  threshold  depends  on  various  parameters 
such  as  gas  pressure  and  composition  as  well  as  on  the  fo¬ 
cusing  conditions.  Furthermore,  the  threshold  intensity 
can  vary  by  as  much  as  two  orders  of  magnitude  due  to 
the  presence  of  aerosols  or  dust  in  the  gas  [4] .  For  clean, 
dust-free  air,  the  breakdown  intensity  threshold  is  of  the 
order  of  1015  W / m2 .  Because  of  such  a  high  optical  inten¬ 
sity  required  for  the  initiation  of  the  breakdown,  break¬ 
down  experiments  are  typically  conducted  using  tightly 
focused  laser  beams.  The  resulting  dense  plasma  plumes, 
under  typical  atmospheric  conditions,  are  less  than  one 
centimeter  long. 

In  the  last  decade,  a  different  mechanism  of  free-charge 
generation  in  gases,  through  femtosecond  laser  filamenta¬ 
tion,  has  been  extensively  studied  [5]— [T] .  In  a  femtosec¬ 
ond  laser  filament,  plasma  generation  via  multi-photon 
ionization  of  the  gas  molecules  constitutes  the  natural 
mechanism  that  balances  the  transverse  collapse  of  the 
beam  due  to  self-focusing. 

Femtosecond  laser  filamentation  allows  for  a  substan¬ 
tial  degree  of  control  over  plasma  generation  in  the  am¬ 
bient  air.  The  longitudinal  position  of  the  plasma  fila¬ 


ment  can  be  straightforwardly  manipulated  by  introduc¬ 
ing  a  negative  temporal  chirp  onto  the  laser  pulse.  As 
the  chirped  pulse  propagates  in  the  dispersive  air,  it  is 
gradually  compressed  in  time  domain,  reaching  its  max¬ 
imum  peak  power  at  a  particular  propagation  distance 
determined  by  the  value  of  the  chirp.  The  placement  of 
the  plasma  channels  within  the  transverse  beam  profile 
can  be  controlled  by  beam  shaping  [8]. 

Longitudinally  extended  plasma  channels  generated  in 
air  through  femtosecond  laser  filamentation  have  numer¬ 
ous  potential  applications  such  as  the  guidance  of  elec¬ 
tricity  and  emission  of  THz  radiation.  However,  the  den¬ 
sity  of  free  electrons  generated  through  femtosecond  fila¬ 
mentation  in  air,  under  typical  conditions,  is  only  about 
1022  m-3  [9],  which  is  too  low  for  many  applications. 

In  this  Letter,  we  report  experiments  on  the  gener¬ 
ation  of  dense  plasma  channels  through  the  combined 
application  of  ultraintense  femtosecond  and  high-energy 
nanosecond  laser  pulses.  Following  the  terminology  used 
in  the  laser  particle  acceleration  community,  we  will 
refer  to  this  approach  as  the  “igniter- heater”  scheme. 
In  this  approach,  the  femtosecond  igniter  pulse  creates 
seed  free  electrons  via  multi-photon  ionization.  The 
co-propagating  nanosecond  heater  pulse  multiplies  these 
seed  electrons  in  an  optically-driven  avalanche.  Our  goal 
is  to  develop  a  viable  technique  for  the  controllable  gener¬ 
ation  of  extended  and  dense  plasma  channels  and  channel 
arrays  in  air,  which  are  of  interest  in  many  applications. 

The  combined  femtosecond-nanosecond  excitation  of 
plasma  in  gases  has  been  previously  used  for  the  gener¬ 
ation  of  short,  dense  plasma  channels  used  in  table-top 
particle  acceleration  [10]  and  in  Raman  amplification  of 
laser  light  in  plasma  channels  [11].  The  experimental  ge¬ 
ometries  used  in  those  works  resulted  in  the  generation 
of,  at  most,  several  centimeter- long  dense  plasma  chan¬ 
nels  and  were  not  straightforwardly  applicable  for  the 
generation  of  such  channels  at  range. 

In  the  context  of  filamentation  science,  experiments  on 
dual-pulse  femtosecond-nanosecond  excitation  have  been 
also  reported  previously.  However,  the  emphasis  has  been 
on  heating,  maintenance  or  revival  of  the  dilute  plasma 
in  filaments,  not  on  the  substantial  multiplication  of  the 
number  of  electrons  in  an  avalanche- like  process  driven 
by  the  nanosecond  pulse  [12], [13].  In  the  above  works, 
the  femtosecond  and  nanosecond  pulses  in  the  two-pulse 
sequence  were  separated  by  a  long  time  interval  of  up  to 
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several  milliseconds.  In  such  a  situation,  the  nanosecond 
pulse  causes  the  liberation  of  the  small  amount  of  elec¬ 
trons  that  attach  themselves  to  neutral  oxygen  molecules 
in  the  air  filament.  The  femtosecond  filament  can  be  re¬ 
vived  in  this  scenario,  but  plasma  density  in  the  revived 
filament  typically  remains  low. 

In  the  present  work,  plasma  generation  occurs  in  a 
different  regime,  in  which  the  nanosecond  heater  pulse 
directly  amplifies  free  electrons  left  in  the  wake  of  the 
femtosecond  igniter  pulse  that  undergoes  filament  at  ion. 
Accordingly,  the  two  pulses  in  our  case  are  temporally 
overlapping  or  nearly  overlapping.  As  we  will  show,  our 
approach  results  in  a  substantial  boost  of  the  plasma  den¬ 
sity  relative  to  that  in  the  femtosecond  filament  alone. 
The  realization  of  our  approach  requires  heater  pulses 
with  high  pulse  energies  of  the  order  of  several  Joules. 
This  requirement  may  have  been  the  reason  why  similar 
experiments  have  not  been  reported  previously. 

Let  us  make  some  basic  estimates.  Complete  single¬ 
electron  ionization  of  oxygen  and  nitrogen  in  the  air  fil¬ 
ament  will  require  certain  energy  that  the  nanosecond 
heater  pulse  will  have  to  expend  per  unit  length  of  the 
generated  plasma  channel.  The  ionization  potentials  of 
oxygen  and  nitrogen  molecules  are  approximately  12  eV 
and  16  eV,  and  their  fractional  contents  in  the  common 
air  are  21%  and  78%,  respectively.  The  total  particle 
density  in  the  air  is  2.7xl025  m-3.  The  diameter  of  a 
typical  femtosecond  laser  filament  in  air  is  about  100  fi m. 
The  ionization  of  all  molecules  in  such  a  channel  will  re¬ 
quire  the  expenditure  of  about  0.5  J  of  energy  per  meter 
of  the  channel.  Contemporary  commercial  Q-switched 
Nd:YAG  lasers  are  available  with  pulse  energies  of  up  to 
50  J.  Thus,  from  the  energy  standpoint,  the  generation 
of  a  completely  ionized,  100  /im-thick  plasma  channel  or 
an  array  of  several  such  channels,  with  the  length  of  the 
order  of  tens  of  meters,  is  feasible. 

To  estimate  the  peak-power  requirement  for  the  heater 
pulse,  we  will  follow  [14]  and  assume  that,  in  the  laser- 
driven  avalanche,  the  energy  of  a  free  electron  in  the  laser 
field  is  increased  by  one  ponderomotive  energy  on  ev¬ 
ery  collision  with  a  neutral  or  already  ionized  molecule. 
Through  the  entire  avalanche  process,  excluding  few  fi¬ 
nal  electron  generations,  the  densities  of  both  electrons 
and  ions  remain  low  compared  to  the  density  of  neutral 
molecules.  Under  these  conditions,  collisions  between 
electrons  and  ions  can  be  ignored  as  they  are  much  less 
frequent  compared  to  the  electron  collisions  with  neutral 
molecules. 

The  rate  equation  for  the  energy  of  a  particular  elec¬ 
tron  reads: 


de(t)  _  e2I(t) 
dt  ~  2 ce0meL02  '  Jem[>  ’ 


(1) 


where  I (t)  is  the  time-dependent  intensity  of  the  nanosec¬ 
ond  heater  pulse,  e  and  me  are  electron  charge  and  mass, 
respectively,  c  is  the  speed  of  light  in  free  space,  eo  is  the 
dielectric  constant  of  vacuum,  and  ujq  is  the  optical  fre¬ 
quency  of  the  heater  pulse.  fem  is  the  frequency  of  colli¬ 


sions  between  electrons  and  neutral  molecules.  It  can  be 
estimated  as  follows: 

~  & erri  ' N m  ’  ^e(^)  =  ^em  -^m  \/ 771  e  ,  (2) 


where  aern  is  the  collisional  cross-section  which  is  as¬ 
sumed  to  be  independent  of  the  electron  energy.  Nm 
is  the  density  of  molecules  in  the  air,  which  remains  ap¬ 
proximately  constant  through  the  entire  development  of 
the  avalanche,  except  for  the  last  few  generations.  ve(t) 
is  the  time-dependent  speed  of  the  electron. 

With  the  exception  of  the  seed  electrons  provided  by 
the  femtosecond  filament,  the  energy  evolution  for  each 
particular  electron,  irrespectively  to  which  electron  gen¬ 
eration  in  the  avalanche  it  belongs  to,  starts  from  a  low 
value  of  the  order  of  one  ponderomotive  energy  and  ends 
with  the  energy  approximately  equal  to  the  ionization  po¬ 
tential  E  of  a  neutral  molecule.  Equation  (1)  is  straight¬ 
forwardly  integrated  from  0  to  r,  which  is  the  time  inter¬ 
val  required  for  an  electron  to  gain  energy  equal  to  the 
ionization  potential.  (We  assume,  for  simplicity,  that  all 
molecules  in  the  air  have  approximately  the  same  ioniza¬ 
tion  potential.)  The  result  can  be  written  as  follows: 


r 


2ce0meujl 


y/2meE, 


(3) 


where  Iq  is  the  peak  intensity  of  the  nanosecond  heater 
pulse,  assumed  to  have  a  rectangular  temporal  profile, 
for  simplicity. 

Through  the  development  of  the  avalanche,  the  total 
number  of  electrons  Ne(t)  approximately  doubles  with 
every  new  generation  of  electrons.  The  doubling  period 
equals  r  given  by  the  formula  (3)  above.  Thus 


Ne(t)  =  N0  •  2t/T  ,  (4) 


where  Nq  is  the  density  of  the  seed  electrons.  At  break¬ 
down,  all  molecules  in  the  air  should  be  ionized  after  the 
passage  of  the  heater  pulse,  i.e.  Ne(T )  «  ATm,  where  T 
is  the  duration  of  the  heater  pulse.  Combining  this  with 
equation  (3),  we  find  the  breakdown  threshold  intensity 
of  the  heater  pulse  in  the  form: 

_  2ce0meuj^y/2m^E  ln(Nm/N0 ) 
th  ~  e2aemNmT  '  ln2  W 

The  ratio  of  logarithms  in  the  above  equation  equals  the 
approximate  number  of  electron  generations  required  to 


500fs  igniter  pulse  Aperture 


Lens  Beam 
combiner 


Plasma 

collection 

electrodes 


Plasma 

channel 


Axicon 


10ns  heater  pulse 


FIG.  1:  Experimental  setup  used  in  the  igniter- heater  excita¬ 
tion  scheme  in  air. 
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FIG.  2:  Left:  Photograph  of  a  dilute  plasma  channel  gener¬ 
ated  through  filamentation  of  the  femtosecond  igniter  pulse 
alone.  To  make  the  faint  filament  visible,  the  image  is  ob¬ 
tained  by  integration  over  100  laser  shots.  Right:  Single-shot 
image  of  a  dense  plasma  channel  generated  through  the  joint 
application  of  the  igniter  and  heater  pulses.  The  inset  shows 
the  fragmentation  of  the  channels  into  discrete  plasma  bub¬ 
bles. 


entirely  ionize  the  air.  For  a  spontaneous  (unseeded) 
optical  breakdown  of  common  air,  it  is  estimated  that 
this  number  is  about  40  [3],  which  corresponds  to  the 
seed  electron  density  of  the  order  of  1013m-3.  Typical 
electron  density  in  a  femtosecond  laser  filament  in  the 
air  is  about  1022  m-3  [9].  If  the  avalanche  is  seeded  by 
a  filament,  the  number  of  electron  generations  required 
for  the  complete  ionization  will  be  about  11.  Thus  even 
though  the  dependence  of  the  breakdown  threshold  in¬ 
tensity  on  the  seed  electron  density  is  very  weak  (logary- 
thmic),  seeding  by  the  filament  will  reduce  the  threshold 
breakdown  intensity  for  the  heater  pulse  by  a  factor  of 
about  4,  a  substantial  improvement  relative  to  the  un¬ 
seeded  breakdown. 

For  the  experimental  demonstration  of  our  approach 
we  use  the  setup  shown  schematically  in  Figure  1.  The 
optical  breakdown  of  ambient  air  is  driven  by  lOns-long 
heater  pulses  generated  by  a  Q-switched  Nd:YAG  laser 
operating  at  1064  nm  wavelength.  The  maximum  pulse 
energy  of  the  heater  laser  is  3.3  J.  The  heater  pulses 
are  combined  with  femtosecond  igniter  pulses  at  800  nm 
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FIG.  3:  Plasma-density  enhancement  factor  relative  to  the 
density  in  the  femtosecond  igniter  filament,  vs.  time  delay 
between  the  igniter  and  heater  pulses.  Both  pulses  are  at 
their  maximum  energies  (15  mJ  and  3.3  J  for  the  igniter  and 
heater  pulses,  respectively). 
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FIG.  4:  Plasma-density  enhancement  factor  as  a  function  of 
the  heater-pulse  energy,  at  the  optimum  delay  between  the 
igniter  and  heater  pulses. 


wavelength.  The  igniter  pulses  are  temporally  chirped 
to  about  500  fs  duration  in  order  to  avoid  filamentation 
and  associated  damage  in  the  dichroic  beam  combiner. 
For  the  purposes  that  will  be  clear  shortly,  the  heater 
beam  is  focused  by  an  axicon  lens  with  the  apex  angle 
of  175°.  The  femtosecond  igniter  is  weakly  focused  by 
an  ordinary  lens  with  the  focal  length  of  0.5  m.  Prior 
to  the  lens,  the  igniter  beam  is  apertured  in  order  to 
clear  a  hole  made  on  the  axis  of  the  axicon  lens  used  to 
focus  the  heater.  The  maximum  energy  of  the  igniter 
pulse  in  the  filamentation  zone  is  about  15  mJ.  Prior  to 
focusing,  the  diameters  of  the  incident  heater  beam  and 
the  apertured  igniter  beam  are  20  mm  and  5  mm,  respec¬ 
tively.  The  pointing  of  the  igniter  beam  is  adjusted  in 
order  to  overlap  the  femtosecond  igniter  filament  with 
the  ^30  cm- long  linear  focal  zone  of  the  axicon.  In  our 
focusing  geometry,  the  waist  of  the  igniter  pulse  is  ap¬ 
proximately  in  the  middle  of  the  linear  focus  zone  for  the 
heater.  Both  igniter  and  heater  lasers  operate  at  10  Hz 
pulse  repetition  frequency.  The  pulsed  outputs  from  the 
two  lasers  are  synchronized  to  better  than  1  ns  accuracy 
through  a  common  triggering.  The  time  delay  between 
the  pulses  can  be  varied  by  means  of  an  electronic  delay 
unit  placed  in  the  trigger  path  for  the  heater  laser. 

With  the  heater  beam  turned  off,  the  femtosecond  ig¬ 
niter  creates  an  ^20  cm- long  dilute  filament,  a  photo¬ 
graph  of  which  is  shown  in  the  left  part  of  Figure  2.  To 
make  this  faint  filament  visible,  the  photograph  was  inte¬ 
grated  over  about  100  laser  shots.  With  both  igniter  and 
heater  pulses  applied  at  their  maximum  available  ener¬ 
gies  and  with  the  temporal  delay  between  the  pulses  op¬ 
timized,  a  seeded  optical  breakdown  is  observed  reliably 
on  every  laser  shot.  The  photograph  of  the  breakdown 
region  is  shown  in  the  right  part  of  Figure  2.  In  this 
case,  the  plasma  channel  is  bright  and  clearly  visible  on 
the  single-shot  image.  The  length  of  the  densified  plasma 
channel  is  about  10  cm.  The  channel  is  fragmented  into 
multiple  plasma  bubbles  separated  by  about  1mm  and 
positioned  irregularly  along  the  channel.  This  fragmen¬ 
tation  is  a  result  of  the  local  plasma  shielding  for  the 
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FIG.  5:  Dependence  of  the  plasma  density  generated  through 
the  joint  application  of  the  igniter  and  heater  pulses,  on  the 
density  of  the  seed  plasma.  Both  seed  and  densified  plasma 
densities  are  plot  on  the  same  arbitrary  unit  scale. 

heater  pulse.  Using  a  heater  beam  in  the  form  of  the  self- 
healing  Bessel  beam  allows  for  the  reconstruction  of  the 
main  intensity  peak  on  the  beam  axis  after  this  peak  is 
locally  scattered  by  the  dense  plasma  generated  through 
the  joint  application  of  the  igniter  and  heater  pulses. 

No  breakdown  is  observed  when  the  heater  beam  is 
applied  at  its  maximum  available  energy  but  the  igniter 
beam  is  turned  off. 

In  order  to  quantify  the  density  of  the  generated 
plasma,  on  a  relative  scale,  we  use  a  simple  capacitive 
plasma  probe  described  in  detail  elsewhere  [15].  In  this 
particular  case,  the  probe  has  1  cm  x  1  cm  collection  elec¬ 
trodes  that  are  separated  by  1.5  cm,  charged  to  100  V, 
and  placed  approximately  in  the  middle  of  the  densified 
plasma  channel. 

Figure  3  shows  the  dependence  of  the  plasma  density 
generated  through  the  joint  application  of  the  igniter  and 
heater  pulses,  as  a  function  of  the  delay  between  the 
pulses.  Positive  values  of  the  delay  correspond  to  the 
nanosecond  heater  pulse  lagging  behind  the  femtosecond 
igniter.  Both  pulses  are  at  their  maximum  available  en¬ 
ergies.  The  plasma  density  is  shown  relative  to  the  den¬ 
sity  with  the  heater  pulse  turned  off  (i.e.  relative  to  the 
density  of  plasma  in  the  seed  filament).  The  maximum 
plasma-enhancement  factor  is  about  200.  For  the  delays 


exceeding  ^10  ns,  the  two  pulses  are  no  longer  tempo¬ 
rally  overlapping,  and  the  breakdown  is  seeded  by  the 
electrons  left  in  the  wake  of  the  femtosecond  pulse. 

The  dependence  of  the  plasma  density  enhancement 
factor  on  the  energy  of  the  heater  pulse  is  shown  in 
Figure  4.  The  energy  of  the  igniter  pulse  is  kept  at  a 
maximum,  and  the  temporal  delay  between  the  pulses 
is  optimized  to  maximize  the  signal.  A  clear  thresh¬ 
old  with  respect  to  the  energy  of  the  heater  pulse  is 
evident,  as  expected.  The  threshold  heater  energy  of 
about  2  J  corresponds  to  the  estimated  peak  on-axis  in¬ 
tensity  in  the  linear  focus  zone  of  the  heater  beam  of 
about  1.5x  1015  W/m1 2 3 4 5 6 7 8. 

In  Figure  5,  we  show  the  dependence  of  the  gener¬ 
ated  electron  density  on  the  density  of  the  seed  electrons 
provided  by  the  igniter  filament.  The  density  of  the  seed 
electrons  was  varied  by  changing  the  energy  of  the  igniter 
pulse.  The  arbitrary  units  in  both  vertical  and  horizontal 
axes  of  this  plot  are  the  same.  The  threshold  with  respect 
to  the  density  of  the  seed  electrons  is  evident.  The  pres¬ 
ence  of  this  threshold  appears  to  be  counter-intuitive. 
It  is  not  immediately  obvious  why  the  reduction  of  the 
density  of  seed  electrons  by,  say,  a  factor  of  2,  completely 
halts  the  breakdown  process.  An  adequate  model  utiliz¬ 
ing  kinetic  equation  formalism  is  necessary  to  understand 
this  behavior. 

In  conclusion,  we  have  conducted  experiments  on  the 
generation  of  dense  plasma  channels  in  ambient  air  us¬ 
ing  the  dual-pulse  femtosecond-nanosecond  excitation 
scheme.  Reliable  triggering  of  the  optical  breakdown 
by  the  femtosecond  seed  filament  was  demonstrated.  A 
maximum  plasma  enhancement  factor  of  about  200  in 
the  dense  plasma  channel  of  about  10  cm  in  length  was 
measured.  Further  development  of  our  approach  and  its 
scaling  to  higher  energies  of  the  heater  pulse  may  yield 
a  viable  technique  for  the  controllable  generation  of  ex¬ 
tended  plasma  channels  and  channel  arrays  in  the  atmo¬ 
sphere  that  can  be  useful  in  many  applications. 
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We  report  experimental  results  on  the  self-focusing  dynamics  and  filamentation  of  ultra-intense  femtosecond 
Bessel  beams  of  various  orders  in  water.  Regular  filament  patterns  are  formed  on  the  bottle-like  dominant 
intensity  features  of  the  beams.  The  peak-power  thresholds  for  filamentation,  at  a  particular  distance,  are 
relatively  accurately  estimated  by  the  adaptation  of  the  Marburger’s  formula  derived  earlier  for  Gaussian 
beams.  The  nonlinear  conversion  of  the  incident  conical  waves  into  the  localized  spatial  wave  packets 
propagating  near  the  beam  axis  is  observed.  ©  2011  Optical  Society  of  America 
OCIS  codes:  190.5940,  190.7110,  050.4865 


Ultrafast  laser  filamentation  is  an  interdisciplinary  field 
of  modern  physics  that  deals  with  the  propagation  of 
ultra- intense  laser  pulses  in  transparent  dielectrics  [1]- 
[3].  Most  of  the  practical  applications  of  filamentation 
are  based  on  the  properties  of  laser  filaments  in  gases. 
However,  transparent  condensed  media  offer  a  conve¬ 
nient  platform  for  experimental  studies  of  filamentation 
and  related  phenomena,  because  the  power  threshold  for 
self-focusing  in  this  case  is  three  orders  of  magnitude 
lower  than  that  in  gases. 

Early  on  in  filamentation  studies,  it  was  realized  that 
beam  shaping  offers  an  effective  means  of  control  over 
the  extent  of  the  generated  filaments  and  their  place¬ 
ment  within  the  beam  profile.  The  application  of  various 
beam  shapes  to  filament  formation,  in  both  gaseous  and 
condensed  media,  has  been  extensively  studied  in  the 
past.  Examples  of  beams  previously  used  for  the  forma¬ 
tion  of  filaments  with  particular  properties  include  fun¬ 
damental  Bessel  beams  [4,5],  vortex  beams  [6,7],  neck¬ 
lace  beams  [8],  and  self-bending  Airy  beams  [9,10]. 

In  this  Letter,  we  report  experiments  on  filamentation 
of  femtosecond  Bessel  beams  of  higher  order  in  water. 
In  addition  to  being  optical  vortices,  these  beams  are 
diffraction-resistant.  They  have  bottle-like  transverse  in¬ 
tensity  distributions  that  retain  their  shapes  on  propa¬ 
gation  over  extended  distances.  The  goal  of  our  study 
is  to  use  these  properties  of  higher-order  Bessel  beams 
for  the  generation  of  extended,  bottle-like  distributions 
of  filaments  that  may  be  useful  in  applications. 

Ideal  higher-order  Bessel  beams  have  their  transverse 
amplitudes  proportional  to  Jn  (27rr/ro)  exp  (in#),  where 
Jn  is  the  n^-order  Bessel  function  of  the  first  kind,  n  is 
the  beam  order,  (r,  0)  are  transverse  polar  coordinates, 
and  ro  is  a  scale  factor.  These  beams  are  generated  from 
an  input  plane  wave  through  the  phase-only  beam  am¬ 
plitude  modulation  with  the  following  transmission  func¬ 
tion  [11]: 

Tn(r,  9)  =  exp  (■ in6 )  exp  ■  (1) 


Ideal  higher-order  Bessel  beams  have  infinite  energy 
flux  and  therefore  they  cannot  be  experimentally  real¬ 
ized.  The  truncated  versions  of  these  beams  are  obtained 
by  the  application  of  the  above  spatial  amplitude  mod¬ 
ulation  to  a  finite-sized  beam  such  as  a  Gaussian  beam. 
The  finite-energy  Bessel  beams  retain  the  dispersion- 
resistance  property  of  the  ideal  Bessel  beams,  but  only 
for  a  limited  propagation  distance  called  the  linear  fo¬ 
cal  zone.  The  extent  of  this  zone  depends  on  the  input 
beam  parameters.  Throughout  this  Letter,  we  will  refer 
to  the  finite-energy  Bessel  beams  as  Bessel  beams,  al¬ 
though  calling  them  Bessel-Gauss  beams  may  be  more 
appropriate.  Examples  of  numerically  simulated  inten¬ 
sity  profiles  for  selected  higher-order  Bessel  beams  are 
shown  in  Figure  1.  Generally,  Bessel  beams  with  an  order 
of  one  or  larger  have  zero  on-axis  intensity.  The  diameter 
of  their  dominant  intensity  ring  increases  with  the  beam 
order. 

The  experimental  setup  used  for  the  studies  of  filamen¬ 
tation  of  higher-order  Bessel  beams  in  water  is  schemat¬ 
ically  shown  in  Figure  2.  The  input  Gaussian  beam  with 
18  mm  beam  diameter  is  generated  by  a  commercial 
ultrafast  Ti: Sapphire-based  amplifier  that  operates  at 
10  Hz  pulse  repetition  frequency  and  outputs  40fs-long 
pulses  with  a  center  wavelength  of  800  nm.  The  input 
beam  is  transformed  into  Bessel  beams  of  various  orders 
through  the  spatial  amplitude  modulation  in  the  form 
(1)  imposed  by  a  computer-controlled,  reflective  Spatial 


Fig.  1.  Transverse  intensity  distributions  for  selected 
higher-order  Bessel  beams  plotted  on  the  same  trans¬ 
verse  scale. 
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Fig.  2.  Schematic  of  the  experimental  setup.  The  inset 
shows  the  word  “Tucson”  that  was  printed  on  the  surface 
of  the  prism  and  imaged  onto  the  camera.  From  this 
image,  the  imaging  depth  is  estimated  at  about  1  mm. 


Light  Modulator  (SLM,  model  X10468  by  Hamamatsu 
Photonics).  The  SLM  matrix  has  600x800  pixels  and 
its  clear  aperture  is  12  mm  by  16  mm,  thus  it  slightly 
clips  the  incident  beam  in  both  horizontal  and  vertical 
dimensions.  The  SLM  did  not  show  any  signs  of  damage 
through  the  continuous  exposure  to  the  incident  laser 
beam  at  the  highest  pulse  energy  level  of  3mJ. 

The  self-focusing  and  filamentation  of  the  Bessel 
beams  takes  place  in  a  glass-walled  cuvette  filled  with 
distilled  water.  The  size  of  the  cuvette  is  1.6  mm  x 
2.0  mm  x  7.6  mm,  and  its  length  is  chosen  to  be  approx¬ 
imately  equal  to  the  dispersion  length  of  the  laser  pulses 
in  water.  The  distance  between  the  entrance  window  of 
the  cuvette  and  the  SLM  is  16  cm.  For  the  particular  re¬ 
alizations  of  higher-order  Bessel  beams  that  we  use,  the 
entire  length  of  the  cuvette  is  well  within  the  linear  focal 
zone  of  the  beams. 

To  image  the  transverse  intensity  distribution  of  the 
beam  propagating  in  the  water,  we  immerse  a  small  glass 
prism  into  the  beam  path.  The  uncoated  glass  surface  of 
the  prism  reflects  about  1.3%  of  the  incident  s-polarized 
beam.  The  reflected  light  is  imaged  onto  a  digital  cam¬ 
era  and  photographed.  For  the  imaging  to  be  linear,  it  is 
important  that  the  reflection  off  the  prism  surface  is  suf¬ 
ficiently  small.  The  prism  is  translated  along  the  beam 
path  which  allows  to  visualize  the  beam  evolution  on 
propagation  through  the  water. 

In  Figure  3,  we  show  the  transverse  intensity  patterns 
imaged  at  5  mm  distance  from  the  entrance  to  the  water 
cell,  for  selected  Bessel  beam  orders.  The  radial  part  of 


Fig.  3.  Selected  higher-order  Bessel  beams  at  low  inten¬ 
sity.  The  ring  intensity  features  of  the  beams  are  modu¬ 
lated  due  to  the  2tt  -  modulo  phase  mask. 


the  phase  modulation  imposed  by  the  SLM,  for  all  cases 
shown,  is  such  that  the  total  phase  shift  from  the  center 
of  the  beam  to  the  edge  equals  100-7T.  The  azimuthal 
component  of  the  modulation  varies  depending  on  the 
beam  order  according  to  (1).  The  peak  power  of  the  laser 
pulses  is  substantially  below  the  self-focusing  threshold, 
thus  the  beams  propagate  in  the  linear  regime. 

It  is  evident  that  the  ring  intensity  features  of  the 
beams,  which  are  supposed  to  be  smooth  for  the  case  of 
ideal  Bessel  beams,  are  modulated.  There  are  exactly  n 
intensity  beads  (or  hot  spots)  along  the  circumference 
of  the  ring  intensity  features  of  the  beam  of  n ^  order. 
The  broadband  nature  of  the  incident  light  coupled  to 
the  wavelength  dependence  of  the  SLM,  the  offset  in  the 
calibration  wavelength  of  the  SLM  with  respect  to  the 
center  wavelength  of  the  laser  source,  as  well  as  a  small 
unmodulated  content  in  the  beam  may  all  be  contribu¬ 
tors  to  this  beeded  beam  structure.  We  have  confirmed, 
however,  that  in  our  case  the  beading  is  predominantly 
caused  by  the  fact  that  the  modulation  phase  imposed  by 
the  SLM  is  defined  modulo  2i r.  The  azimuthal  component 
of  the  phase  modulation  for  the  n^-order  beam  experi¬ 
ences  exactly  n  27r-jumps  along  the  circumference  of  the 
beam.  Thus  the  phase  is  not  defined  at  these  n  equidis- 
tantly  spaced  jump  points,  causing  the  resulting  beam 
pattern  to  have  n  regularly  spaced  hot  spots  on  its  dom¬ 
inant  intensity  features.  The  peak-to-peak  amplitude  of 
the  intensity  modulation  of  the  dominant  ring,  relative 
to  its  mean  value,  is  between  10%  and  20%,  depending 
on  the  beam  order. 

When  the  pulse  energy  is  increased,  distinct  filaments 
start  forming  on  the  dominant  intensity  rings  of  the 
beams.  The  self- focusing  collapse  is  accompanied  by  the 
emission  of  the  forward-propagating  white-light  contin¬ 
uum  radiation  [1]-  [3].  In  the  top  row  of  Figure  4,  the 
white-light  continuum  is  imaged  at  the  pulse-energy  lev¬ 
els  when  filamentation  starts  to  occur  at  a  5  mm  distance 


Fig.  4.  Top  row:  Images  of  the  white  light  emission  at 
the  pulse-energy  levels  corresponding  to  the  onset  of  fil¬ 
amentation  at  5  mm  from  the  entrance  to  the  water  cu¬ 
vette,  for  selected  higher-order  Bessel  beams.  Bottom 
row:  Same  at  the  energy  levels  corresponding  to  the  onset 
of  filamentation  in  the  secondary  ring  intensity  features. 
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Fig.  5.  Critical  power  for  filamentation  for  higher-order 
Bessel  beams.  Experimental  data  points  are  shown  with 
circles.  The  solid  line  is  the  result  of  the  calculation  based 
on  the  adaptation  of  the  Marburger’s  formula. 


from  the  entrance  to  the  water  cell.  These  images  are 
taken  through  a  color-glass  filter  that  allows  visible  light 
to  pass  through  but  blocks  the  pump  light  at  800  nm. 

If  the  dominant  rings  of  the  Bessel  beams  were  per¬ 
fectly  smooth,  2n  +  1  filaments  would  be  expected  to 

th 

form  on  the  main  ring  of  the  n  -order  beam  due  to  the 
azimuthal  modulation  instability  [7,12].  However,  in  our 
case,  the  placement  of  filaments  is  forced  by  the  beaded 
structure  of  the  beams,  causing  the  n^-order  beam  to 
have  exactly  n  filaments  on  its  dominant  ring.  As  the 
pulse  energy  is  increased  further,  filaments  start  forming 
on  the  secondary  rings  of  the  beams,  as  shown  in  the 
bottom  row  of  Figure  4. 

For  bell-shaped  beams,  the  peak  power  threshold  Pth 
for  filamentation,  at  a  particular  distance  Lc,  can  be  rel¬ 
atively  accurately  estimated  by  using  the  semi-empirical 
Marburger’s  formula  [13]: 

T  _ Q.36  77rn0Wo _ 

AoVKPtft/Pcr)1/2  -  0.852]2  -  0.0219  ’ 

where  wq  is  the  input  beam  radius,  Ao  is  the  wavelength 
of  the  pump  beam  in  vacuum,  no  is  the  linear  refractive 
index  of  the  medium,  and  Pcr  is  the  critical  power  for 
self-focusing  at  infinity.  This  formula  can  be  adapted  for 
the  estimation  of  the  filamentation  thresholds  for  higher- 
order  Bessel  beams  at  a  particular  propagation  distance. 
In  our  case,  the  power  in  the  dominant  intensity  ring 
of  the  n^-order  beam  is  approximately  equally  parti¬ 
tioned  between  the  n  formed  filaments,  and  wo  approxi¬ 
mately  equals  to  the  half-width  of  the  ring.  The  fraction 
of  power  in  the  center  ring  is  deduced  from  the  numeri¬ 
cal  solution  of  the  Fresnel  diffraction  integral.  The  result 
of  the  calculation  following  this  procedure,  for  the  case 
of  Lc  =  5  mm,  is  shown  in  Figure  5  with  a  solid  line. 
The  experimental  data  for  selected  beam  orders  is  in  a 
reasonably  good  agreement  with  the  calculation. 

As  the  self-focusing  higher-order  Bessel  beam  propa¬ 
gates  further  into  the  water  cell,  the  intensity  ring  fea¬ 
tures  of  the  beam  self-focus  as  a  whole,  at  the  same  time 


Fig.  6.  Jg  (top  row)  and  J12  (bottom  row)  beams  im¬ 
aged  at  50  mm  from  the  entrance  to  the  water  cuvette, 
at  various  pulse  energies. 


as  the  breakup  of  the  rings  into  individual  filaments  oc¬ 
curs.  In  addition,  the  nonlinear  conversion  of  the  con¬ 
ical  waves  from  the  incident  beam  into  various  spatial 
wavepackets  propagating  near  the  beam  axis  takes  place. 
A  similar  effect  has  been  previously  reported  for  the  case 
of  self- focusing  of  the  fundamental  Bessel  beam  [14].  Ex¬ 
amples  of  these  complex  dynamics  are  shown  in  Figure  6, 
for  two  selected  beam  orders. 

In  summary,  we  have  studied  the  filamentation  of  fem¬ 
tosecond  higher-order  Bessel  beams  in  water.  Regular 
filament  patterns  have  been  found  to  form  on  the  domi¬ 
nant  intensity  rings  of  the  beams.  The  placement  of  fila¬ 
ments  within  the  transverse  beam  profile  was  dominated 
by  the  azimuthal  amplitude  modulation  imposed  by  the 
‘I'K-modulo  phase  mask  used  to  generate  the  beams.  The 
pulse-energy  thresholds  for  the  onset  of  filamentation, 
at  a  particular  propagation  distance,  have  been  rela¬ 
tively  accurately  estimated  by  the  adaptation  of  the 
Marburger’s  formula.  The  nonlinear  conversion  of  the 
incident  conical  waves  into  the  localized  wave  packets 
propagating  near  the  beam  axis  has  been  observed. 

The  authors  acknowledge  support  from  The  United 
States  Air  Force  Office  of  Scientific  Research  under  pro¬ 
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We  report  experiments  and  numerical  simulations  on  supercontinuum  generation  with  femtosecond 
Airy  pulses  in  a  highly  nonlinear  optical  fiber.  The  ability  of  the  Airy  waveform  to  regenerate  its 
dominant  intensity  peak  results  in  the  generation  of  distinct  spectral  features.  Airy  pulses  and  other 
self-healing  temporal  waveforms  may  be  useful  for  the  generation  of  spectra  with  desired  properties. 

PACS  numbers:  42. 65. Re,  42. 81. Dp 


The  dramatic  spectral  broadening  experienced  by  in¬ 
tense  optical  pulses  in  nonlinear  media,  commonly  known 
as  supercontinuum  generation,  is  a  result  of  the  com¬ 
plex  interplay  between  linear  dispersion  and  instanta¬ 
neous  and  delayed  nonlinearity  in  the  medium  [1].  Mi- 
crostructured  optical  fibers  are  particularly  well  suited 
for  supercontinuum  generation  because  the  nonlinearity 
and  dispersion  landscape  in  microstructured  fibers  can 
be  altered  by  the  design  of  the  waveguide  structure  [2] . 

Although  supercontinuum  generation  in  optical  fibers 
has  been  extensively  studied  under  a  vast  variety  of  con¬ 
ditions  in  the  past  (for  a  review,  see,  for  example,  [3]), 
the  overwhelming  majority  of  prior  studies  utilized  in¬ 
tense  optical  pulses  with  bell-shaped  temporal  profiles. 
There  have  been  several  reports  on  the  optimization  of 
supercontinuum  generation  via  pulse  shaping  [4]  -  [6]. 
However,  this  approach  remains  insufficiently  explored. 

Recently,  a  new  kind  of  optical  pulse  has  been  intro¬ 
duced  to  optics  [7].  These  so-called  Airy  pulses  have 
their  field  envelopes  described  in  terms  of  truncated  or 
apodized  Airy  functions  of  time.  These  pulses  are  tempo¬ 
ral  analogues  of  Airy  beams,  also  introduced  in  [7].  Like 
Airy  beams  resist  diffraction  and  self-bend,  Airy  pulses 
resist  dispersion  and  their  dominant  intensity  peaks  ac¬ 
celerate  on  propagation.  Another  important  property  of 
Airy  pulses  is  their  ability  to  regenerate  their  dominant 
intensity  features  should  those  features  be  selectively  at¬ 
tenuated  or  distorted.  All  of  the  above  features  are  linear 
effects.  The  dispersion  resistance  and  self-regeneration 
of  the  temporal  Airy  wavepacket  are  illustrated  by  sim¬ 
ulations  shown  in  Figure  1.  Temporal  Airy  pulses  have 
been  experimentally  realized  and  their  properties  have 
been  studied  in  both  linear  [8]  and  nonlinear  [9]  regimes. 

The  ability  of  Airy  pulses  to  self-heal  makes  them  po¬ 
tentially  useful  in  the  context  of  supercontinuum  genera¬ 
tion.  With  an  input  bell-shaped  pulse  and  in  the  anoma¬ 
lous  dispersion  regime,  which  is  commonly  used  for  the 
generation  of  broadest  continua,  the  spectral  evolution  of 
the  pulse  largely  ends  when  the  pulse  is  transformed  into 
one  or  several  fundamental  solitons.  The  excess  energy  in 
this  case  is  shed  in  the  form  of  dispersive  waves  [10].  The 
placement  of  these  features  in  the  spectrum  is  determined 
by  the  dispersion  landscape  of  the  fiber.  In  the  case  of 
the  self-healing  Airy  pulse,  the  generation  of  new  spec¬ 
tral  components  may  continue  after  the  dominant  peak  of 
the  waveform  has  produced  solitons  and  dispersive  waves 
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FIG.  1:  Top  row:  Input  and  output  temporal  intensity  profiles 
of  a  Gaussian  pulse  propagating  in  a  linear  dispersive  optical 
fiber  that  has  a  second-order  dispersion  only.  The  24  fs- long 
input  pulse  is  severely  dispersed  on  propagation.  Middle  row: 
The  input  pulse  is  an  exponentially  apodized  Airy  pulse  with 
the  same  optical  bandwidth  as  that  of  the  Gaussian  pulse 
above.  Propagation  through  the  same  fiber  leaves  the  domi¬ 
nant  peak  of  the  waveform  largely  intact.  Bottom  row:  The 
dominant  intensity  peak  of  the  input  Airy  pulse  is  cutoff  at 
a  level  of  one-sixth  of  the  peak  value.  The  feature  is  restored 
on  propagation  through  the  same  fiber  as  above. 


and  removed  itself  from  the  continuum-generation  pro¬ 
cess.  Under  the  right  conditions,  the  dominant  feature  of 
the  Airy  waveform  will  be  regenerated.  At  that  point,  it 
will  propagate  with  a  velocity  different  from  that  of  the 
dominant  peak  of  the  original  pulse,  affecting  the  place¬ 
ment  of  the  newly  generated  spectral  components.  Thus 
the  use  of  Airy  pulses  and  other  sophisticated  self-healing 
temporal  waveforms  for  continuum  generation  may  allow 
for  control  over  the  generated  spectrum  by  linear  optic 
means. 

In  this  Letter,  we  report  the  results  of  a  numerical  and 
experimental  study  of  supercontinuum  generation  with 
Airy  pulses  in  a  microstructured  optical  fiber.  Our  re¬ 
sults  show  that  the  self-healing  property  of  the  pulse,  un¬ 
der  certain  conditions,  results  in  nontrivial  cyclic  behav¬ 
ior.  The  dominant  feature  of  the  pulse  creates  a  soliton- 
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FIG.  2:  Top  row:  Input  and  output  temporal  waveforms  for 
the  case  of  the  propagation  of  an  Airy  pulse  with  -60,000  fs3 
cubic  spectral  phase  in  a  1  m-long  nonlinear  fiber  with  the 
parameters  described  in  the  text.  Bottom  row:  Correspond¬ 
ing  spectra.  Solitons  and  dispersive  waves  are  marked  with 
letters  S*  and  D*,  where  the  subscript  indicates  the  order  of 
appearance  of  these  features. 


dispersive  wave  pair,  regenerates  itself  through  the  self- 
healing  action  of  the  Airy  pulse,  then  creates  another 
soliton-dispersive  wave  pair,  and  so  on.  Such  propagation 
dynamics  results  in  the  generation  of  distinct  features  in 
the  output  optical  spectrum. 

An  Airy  pulse  is  created  by  imposing  a  cubic  phase 
onto  the  spectrum  of  an  input  Gaussian  pulse.  The  cu¬ 
bic  phase  can  be  either  positive  or  negative,  which  cor¬ 
responds  to  the  Airy  pulse  with  the  tail  lagging  behind 
the  dominant  peak  of  the  pulse  or  vice  versa.  The  Airy 
pulse  shown  in  Figure  1  has  a  negative  value  of  the  cubic 
spectral  phase.  Accordingly,  the  tail  of  the  waveform  in 
this  case  is  propagating  in  front  of  the  dominant  peak. 

The  orientation  of  the  tail  of  the  pulse  with  respect 
to  the  pulse  waveform  has  two  important  consequences. 
First,  linear  third-order  dispersion,  which  has  a  partic¬ 
ular  sign,  will  add  to  or  subtract  from  the  cubic  phase 
of  the  input  Airy  pulse,  resulting  in  the  extension  or  re¬ 
duction  of  the  dispersion-resistant  propagation  range  of 
the  pulse.  In  fact,  the  generation  of  temporal  Airy  wave¬ 
forms  from  input  Gaussian  pulses,  as  a  result  of  prop¬ 
agation  through  a  fiber  with  third-order  dispersion,  has 
been  described  as  early  as  1979  [11],  much  before  the  con¬ 
cept  of  the  dispersion-resistant  Airy  pulse  had  been  pro¬ 
posed.  Second,  solitons  produced  through  the  collapse 
of  the  dominant  peak  of  the  Airy  pulse  generally  tend  to 
lag  behind  the  main  pulse  waveform.  Thus  for  the  case 
when  the  tail  propagates  ahead  of  the  main  peak  of  the 
pulse,  the  generated  solitons  cleanly  detach  themselves 
from  the  waveform  and  leave  the  self-healing  of  the  pulse 
unaffected.  When  the  tail  of  the  pulse  propagates  be¬ 
hind  its  main  peak,  the  generated  solitons  pass  through 
the  tail  and  interfere  with  the  self-healing. 

The  parameters  of  the  optical  fiber  used  in  the  sim¬ 
ulations  throughout  this  report  closely  match  those  of 
the  fiber  used  in  our  experiments.  The  model  fiber 
has  the  nonlinear  coefficient  7  equal  to  O.ll(W-m)-1, 
zero-dispersion  wavelength  of  745  nm,  and  the  maximum 
group- velocity  dispersion  of  240ps/nm-km,  which  occurs 


at  1.96  /im  wavelength.  The  fiber  dispersion  is  anoma¬ 
lous  for  all  wavelengths  above  745  nm.  At  800  nm,  which 
is  the  center  wavelength  of  the  input  pulses,  the  value  of 
dispersion  is  35  ps/nm-km.  The  dispersion  profile  of  the 
model  fiber  is  fitted  by  a  10-order  polynomial  in  optical 
frequency. 

The  input  Airy  pulses  are  numerically  generated  by  im¬ 
posing  a  cubic  spectral  phase  of  minus  or  plus  60,000  fs3, 
onto  the  spectrum  of  a  24  fs- long,  transform-limited 
Gaussian  pulse.  Under  these  conditions,  the  dominant 
peak  of  the  pulse  contains  about  one  half  of  the  total  en¬ 
ergy  of  the  waveform.  The  pulse  evolution  is  simulated 
by  using  the  numerical  model  adapted  from  [12].  In  the 
model,  the  optical  field  is  assumed  to  remain  linearly  po¬ 
larized  throughout  the  entire  propagation  in  the  fiber. 

Figure  2  shows  the  temporal  and  spectral  outputs  from 
a  1  m-long  nonlinear  fiber  with  the  parameters  listed 
above,  for  the  case  of  an  input  Airy  pulse  with  peak  power 
of  1.5  kW  and  a  negative  cubic  spectral  phase  equal  to 
-60,000 fs3.  The  tail  of  the  waveform  in  this  case  propa¬ 
gates  in  front  of  its  main  peak.  Three  soliton-formation 
events  are  evident  from  the  side-by-side  examination  of 
the  temporal  and  spectral  shapes  of  the  pulse  as  it  prop¬ 
agates  through  the  fiber.  All  three  solitons  are  pro¬ 
duced  through  the  reshaping  of  the  main  intensity  peak 
of  the  Airy  pulse,  which  recreates  itself  after  each  soliton- 
formation  event.  The  spectral  features  corresponding  to 
the  generated  solitons  and  dispersive  waves  are  marked 
in  the  Figure  by  letters  S1-S3  and  D1-D3,  where  the  sub¬ 
script  indicates  the  order  of  appearance  of  these  features. 
After  their  formation,  solitons  and  dispersive  waves  re¬ 
main  essentially  invariant  on  propagation,  except  for  the 
continuous  frequency  red  shifting  of  the  soliton  spectral 
peaks  caused  by  higher-order  dispersion  and  delayed  Ra¬ 
man  nonlinearity. 

The  red-shifted  solitons  gradually  slow  down  and  may 
catch  up  with  dispersive  waves  in  time  domain.  In  that 
case,  the  cross-phase  modulation  induced  by  the  intense 
soliton  field  causes  the  redistribution  of  the  energy  of 
the  dispersive  wave  and  results  in  the  appearance  of  an 
additional  blue-shifted  spectral  peak,  as  described  in  [13]. 
The  left-most  peak  at  570  nm  in  the  bottom  right  panel 
of  Figure  2  is  the  result  of  such  an  interaction  between 
the  first  emitted  soliton  Si  and  the  dispersive  wave  Di. 

The  placement  of  the  spectral  peak  corresponding  to 
a  particular  dispersive  wave  is  determined  by  the  wave- 
number  matching  between  the  dispersive  wave  and  the 
higher-order  soliton  that  undergoes  fission,  in  the  imme¬ 
diate  vicinity  of  the  point  where  the  peak  power  of  the 
soliton  in  the  highest  [14].  For  the  case  of  a  fundamental 
soliton,  the  matching  condition  has  the  following  form: 

n>  2 

where  /3n  (c os)  are  the  higher-order  dispersion  coefficients 
of  the  fiber  at  00  s,  the  optical  frequency  of  the  soliton, 
oj d  is  the  center  frequency  of  the  dispersive  wave,  7  is  the 
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nonlinear  parameter  of  the  fiber,  and  Ps  is  the  soliton 
peak  power.  Although  (1)  is  only  valid  for  fundamental 
solitons,  it  is  frequently  used  for  estimations  of  the  dis¬ 
persive  wave  placement  in  the  cases  of  solitons  of  higher 
order.  The  center  wavelengths  of  the  three  generated  dis¬ 
persive  waves  in  the  case  shown  in  Figure  2  are,  in  order 
of  their  appearance,  600  nm,  668  nm,  and  690  nm.  The 
estimates  based  on  the  formula  (1)  are  580  nm,  639  nm, 
and  663  nm.  These  estimates  are  consistently  lower  than 
the  values  derived  from  the  simulations.  The  poor  agree¬ 
ment  may  be  caused  by,  among  other  things,  the  accel¬ 
erating  nature  of  the  Airy  waveform,  for  which  case  the 
condition  (1)  is  overly  simplified. 

Figure  3  provides  a  closer  look  at  the  first  two  soliton- 
formation  events.  The  top  left  panel  shows  the  temporal 
profile  of  the  pulse  after  it  propagated  6  cm  through  the 
fiber.  At  that  point,  the  main  peak  of  the  input  pulse  has 
transformed  itself  into  a  higher-order  soliton  that  sub¬ 
sequently  underwent  fission  and  shed  a  dispersive  wave. 
The  main  peak  of  the  waveform  is  nearly  regenerated  and 
ready  to  produce  another  soliton.  The  top  right  panel 
shows  the  pulse  profile  after  16  cm  of  propagation.  The 
regenerated  main  peak  of  the  Airy  pulse  has  just  created 
another  soliton. 

In  the  spectral  domain,  the  soliton  and  dispersive 
wave  created  through  the  first  collapse  of  the  main 
Airy  peak  are  straightforwardly  identifiable  at  both  6  cm 
and  16  cm.  Spectral  features  resulting  from  the  second 
soliton-formation  event  are  not  as  obvious.  To  properly 
identify  them,  we  computed  a  Fourier  transform  of  the 
portion  of  the  temporal  waveform  at  16  cm  that  contained 
the  peak  marked  by  S2. 

In  the  above  example,  the  Airy  pulse  had  a  negative 
value  of  the  cubic  spectral  phase,  in  which  case  the  tail 
of  the  Airy  waveform  propagated  in  front  of  its  main 
peak.  The  solitons  and  dispersive  waves  generated  from 
the  main  peak  of  the  pulse  lagged  behind  the  waveform 
and  did  not  interfere  with  its  ability  to  self-regenerate. 
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FIG.  3:  Top  row,  left:  Temporal  waveforms  of  the  Airy  pulse 
with  -60,000  fs3  cubic  spectral  phase  after  6  cm  propagation  in 
the  fiber.  The  main  peak  of  the  Airy  waveform  has  produced  a 
soliton-dispersive  wave  pair  and  is  nearly  regenerated.  Right: 
Same  after  16  cm  propagation.  The  regenerated  main  peak  of 
the  pulse  has  produced  another  soliton-dispersive  wave  pair. 
Bottom  row:  Corresponding  spectra. 
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FIG.  4:  Input  Airy  pulse  with  +60,000  fs3  cubic  spectral  phase 
propagates  in  1  m-long  nonlinear  fiber.  Top  row:  Temporal 
input  and  output  waveforms.  Bottom  row:  Corresponding 
spectra.  Only  one  distinct  dispersive-wave  peak  is  produced. 


The  case  corresponding  to  the  same  value  but  oppo¬ 
site  sign  of  the  cubic  spectral  phase  is  shown  in  Figure  4. 
As  in  the  previous  case,  the  dominant  intensity  peak  of 
the  input  Airy  pulse  relatively  quickly  reshapes  itself  into 
a  higher-order  soliton  which  subsequently  morphs  into  a 
fundamental  soliton  and  emits  a  dispersive  wave.  How¬ 
ever,  in  this  case,  the  generated  soliton  passes  through 
the  tail  of  the  Airy  pulse  and  perturbs  the  phase  of  this 
waveform  via  cross-phase  modulation.  As  a  result  of  that 
perturbation,  the  Airy  waveform  loses  its  ability  to  re¬ 
generate  its  main  peak.  Two  additional  solitons  are  still 
produced  by  the  remaining  part  of  the  pulse,  but  these 
features  originate  from  the  secondary  peaks  of  the  wave¬ 
form,  and  not  from  the  rebuilt  dominant  intensity  fea¬ 
ture.  Only  one  distinct  dispersive-wave  peak  is  evident 
in  the  output  spectrum  in  this  case.  As  in  the  case  with 
negative  cubic  phase,  the  leftmost  feature  in  the  output 
spectrum  results  from  the  interaction  of  the  first  emitted 
soliton  with  its  own  dispersive  wave. 

In  order  to  validate  the  pulse  propagation  scenario 
put  forth  above  by  numerical  simulations,  we  conducted 
experiments  on  supercontinuum  generation  with  Airy 
pulses  in  a  nonlinear  fiber.  In  the  experiments,  the  in¬ 
put  Gaussian  pulses  are  generated  by  a  commercial  mod- 
elocked  Ti: Sapphire  laser  oscillator  operating  at  800  nm 
center  wavelength  and  80  MHz  pulse  repetition  frequency. 
Pulses  out  of  the  oscillator  are  de-chirped  to  about  24  fs 
duration,  which  is  close  to  the  transform  limit,  using  a 
dispersive  prism  pair.  Airy  pulses  are  produced  by  impos¬ 
ing  a  cubic  spectral  phase  onto  these  transform-limited 

Gaussian  pulses,  using  a  commercial  computer-controlled 
'T'A/r 

pulse  shaper  (Silhouette liV1  by  Coherent).  The  maxi¬ 
mum  cubic  phase  attainable  in  our  setup  is  limited  to 
plus  or  minus  60,000 fs3,  which  matches  the  magnitude 
of  the  cubic  phase  used  in  the  simulations.  The  temporal 
profile  of  the  shaped  pulses  is  analyzed  using  a  commer¬ 
cial  FROG  system. 

The  shaped  pulses  are  coupled  into  a  1  m-long  nonlin¬ 
ear  fiber  with  the  parameters  that  closely  match  those 
used  in  the  computer  simulations  discussed  above.  The 
mode-field  diameter  of  the  fiber  is  about  1.4 /im,  its  es- 
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FIG.  5:  Supercontinuum  spectra  recorded  at  the  output  of 
a  1  m-long  microstructured  fiber  with  the  parameters  closely 
matching  those  used  in  the  simulations.  The  case  with  the 
input  Airy  pulse  with  a  negative  cubic  spectral  phase  equal 
-60,000  fs3  (left)  results  in  the  generation  of  three  distinct 
dispersive- wave  features.  The  case  with  the  same  but  op¬ 
posite  spectral  phase  (right)  results  in  the  generation  of  only 
one  distinct  dispersive  wave. 

timated  nonlinear  coefficient  7  equals  105  (W-m)_1,  and 
the  zero  dispersion  wavelength  of  the  fiber  is  745  nm.  The 
coupling  efficiency  of  the  shaped  laser  pulses  into  the  fiber 
in  our  setup  is  about  20%,  which  corresponds  to  0.25  nJ 
of  maximum  pulse  energy  coupled  into  the  fiber.  The 
input  pulse  energy  is  varied  by  adjusting  the  focusing 
of  the  fiber- coupling  arrangement.  The  supercontinuum 
spectra  at  the  output  from  the  fiber  are  recorded  using 
an  optical  spectrum  analyzer. 

Two  examples  of  experimentally  measured  spectra  are 
shown  in  Figure  5.  In  both  cases  shown,  the  energy  of 
input  Airy  pulses  equals  0.12nJ.  The  amount  of  the  cu¬ 
bic  spectral  phase  imposed  by  the  pulse  shaper  in  or¬ 
der  to  generate  these  waveforms  is  -60,000  fs1 2 3 4 5 6 7  (left)  and 
+60,000 fs3  (right).  Under  these  conditions,  the  esti¬ 
mated  peak  power  of  the  Airy  pulses  coupled  into  the 
fiber  is  1.5  kW,  matching  the  value  used  in  the  simula¬ 
tions. 


For  the  case  with  a  negative  cubic  spectral  phase,  the 
generated  supercontinuum  spectrum  has  three  distinct 
and  approximately  equidistant  dispersive  wave  peaks,  in 
qualitative  agreement  with  the  simulation  shown  in  Fig¬ 
ure  2.  For  the  case  with  a  positive  cubic  phase,  only 
one  distinct  dispersive-wave  feature  is  produced,  again 
in  agreement  with  the  simulations.  In  both  cases,  the 
left-most  peak  in  the  spectrum  is  due  to  the  interaction 
between  the  first  emitted  soliton  and  its  own  dispersive 
wave.  The  quantitative  discrepancy  between  the  experi¬ 
mental  results  and  the  simulations  is  likely  due  to  the  in¬ 
complete  knowledge  of  the  fiber  parameters  as  well  as  to 
the  residual  fiber  birefringence,  which  was  not  accounted 
for  in  our  model. 

In  conclusion,  we  conducted  experiments  and  numeri¬ 
cal  simulations  on  supercontinuum  generation  with  self- 
healing  Airy  pulses  in  a  highly  nonlinear  fiber.  The  pulse 
propagation  dynamic  is  found  to  be  different  for  the  cases 
when  the  tail  of  the  Airy  waveform  propagates  in  front 
of  and  behind  the  main  peak  of  the  pulse.  In  the  first 
case,  the  dominant  peak  of  the  pulse  recreates  itself  after 
several  soliton  formation  events,  resulting  in  the  gener¬ 
ation  of  distinct  dispersive-wave  features  in  the  output 
spectrum.  In  the  second  case,  the  soliton  generated  from 
the  first  collapse  of  the  main  peak  of  the  pulse  passes 
through  the  tail  of  the  waveform  and  disrupts  its  abil¬ 
ity  to  self-regenerate.  The  results  of  the  modeling  are  in 
qualitative  agreement  with  the  experiments. 
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We  report  experiments  on  the  generation  of  third  and  fifth  harmonics  of  millijoule-level,  tightly  fo¬ 
cused,  femtosecond  laser  pulses  at  2.2  gm  wavelength  in  air.  The  measured  ratio  of  yields  of  the  third 
and  fifth  harmonics  in  our  setup  is  about  2  •  10-4.  This  result  contradicts  the  recent  suggestion  that  the 
Kerr  effect  in  air  saturates  and  changes  sign  in  ultra-intense  optical  fields.  ©  2011  Optical  Society  of  America 
OCIS  codes:  190.2620,  190.3270,  320.7110 


Linear  and  nearly  instantaneous  dependence  of  the  re¬ 
fractive  index  of  transparent  media  on  the  intensity  of 
the  optical  field  is  the  essence  of  the  electronic  Kerr  ef¬ 
fect,  which  is  the  cornerstone  of  nonlinear  optics.  In  a 
recent  publication,  it  has  been  suggested  that  the  notion 
of  the  instantaneous  electronic  Kerr  effect  must  be  gen¬ 
eralized  [1].  Measurements  of  the  transient  birefringence 
induced  in  various  gases  by  ultraintense  and  ultrashort 
optical  pulses  have  been  interpreted  as  the  saturation 
and  sign  reversal  of  the  nonlinear  refractive  index  at  the 
optical  intensity  levels  on  the  order  of  several  tens  of 
terawatts  per  square  centimeter.  The  common  linear  de¬ 
pendence  of  the  index  of  refraction  on  the  intensity  of 
the  optical  field  has  been  amended  by  the  inclusion  of 
new,  nontrivially  large  terms  proportional  to  second  and 
higher  powers  of  intensity. 

The  field  of  femtosecond  laser  filamentation  is  proba¬ 
bly  the  one  that  would  be  affected  the  most,  should  the 
above  suggestion  be  proven  valid.  For  many  years,  it  was 
believed  that  the  intrinsically  unstable  balance  between 
diffraction  and  self-focusing  in  ultraintense  laser  beams 
undergoing  filamentation  in  gaseous  media  is  stabilized 
by  the  weak  defocusing  action  of  plasma  generated  on 
the  beam  axis  via  multi-photon  ionization  [2]  -  [4].  The 
higher-order  Kerr  terms,  if  they  did  exist,  could  provide 
for  an  alternative  stabilization  mechanism.  The  filament 
propagation  would  be  much  longer  ranged,  as  the  losses 
to  ionization  would  be  reduced  compared  to  the  estab¬ 
lished  filamentation  scenario  based  on  plasma  defocus- 
ing  [5]. 

Effects  associated  with  weak  saturation  of  the  elec¬ 
tronic  Kerr  response  have  been  considered  in  the  con¬ 
text  of  laser  filamentation  previously  [6]  -  [8].  However, 
the  suggestion  that  the  Kerr  effect  may  reverse  sign  and 
cause,  among  other  things,  beam  self-defocusing,  were 
too  radical  a  departure  from  the  established  framework 
of  laser  filamentation  and  nonlinear  optics  in  general.  Al¬ 
most  immediately  after  the  publication  of  [1],  the  results 
reported  in  that  paper  became  the  subject  of  intense  de¬ 
bates. 

Following  [9],  we  will  refer  to  the  claim  put  forth  in  [1] 


as  the  Higher- Order  Kerr  Effect  (HOKE)  theory.  Several 
independent  experimental  tests  of  HOKE  in  the  context 
of  laser  filamentation  have  unambiguously  proven  that 
even  if  the  higher-order  Kerr  terms  do  exist,  they  are  not 
operative  in  common  laser  filaments  in  gases  [10]  -  [12]. 
The  question  still  remains  whether  this  high-order  Kerr 
effect  exists  at  all.  It  has  been  argued  that  the  sign  re¬ 
versal  of  the  Kerr  effect  may  occur  in  a  transient  regime, 
while  neutral  molecules  in  the  gas  are  undergoing  multi¬ 
photon  ionization  [13],  [14]. 

Very  recent  experiments  on  direct  profiling  of  the 
cross-phase  modulation  response  in  a  thin  gas  target  il¬ 
luminated  by  an  ultraintense  and  ultrashort  laser  pulse 
revealed  no  evidence  of  even  a  transient  saturation  of  the 
Kerr  effect,  in  complete  contradiction  with  HOKE  [15]. 
An  alternative  test  that  would  also  look  at  the  pres¬ 
ence  of  a  potentially  transient  HOKE  response  has  been 
proposed  in  [16].  In  this  Letter,  we  report  experimen¬ 
tal  results  on  the  realization  of  that  proposal.  In  agree¬ 
ment  with  [15],  we  found  no  evidence  of  the  higher-order 
Kerr  response.  We  believe  that  at  this  time,  the  com¬ 
bined  body  of  experimental  evidence  is  sufficient  to  per¬ 
manently  disqualify  the  HOKE  paradigm. 

The  experimental  test  reported  here  is  based  on  the 
dramatic  effect  that  the  inclusion  of  the  HOKE  terms 
would  have  on  the  efficiency  of  third  and  fifth  harmonic 
generation  by  a  femtosecond  laser  pulse  with  peak  power 
in  the  vicinity  of  the  alleged  turnover  point  for  the 
higher-order  Kerr  effect.  As  has  been  numerically  shown 
in  [16]  for  tightly  focused  femtosecond  pulses  at  1.3 /im 
wavelength,  the  ratio  of  the  yields  of  the  fifth  and  third 
harmonics  generated  by  such  pulses  is  a  growing  func¬ 
tion  of  the  intensity  of  the  pulse.  This  ratio  saturates  at 
a  value  on  the  order  of  one  if  the  HOKE  terms  are  in¬ 
cluded  in  the  model,  and  it  remains  very  small  10-4) 
when  these  terms  are  omitted.  Such  a  dramatic  differ¬ 
ence  between  the  predictions  of  the  two  models  is  the 
consequence  of  the  fact  that  in  the  established  approach, 
the  fifth  harmonic  generation  is  a  cascade  process  involv¬ 
ing  the  nonlinear  mixing  between  the  already  generated 
third  harmonic  and  the  leftover  pump,  while  in  the  new 
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theory,  the  fifth  harmonic  is  generated  directly  from  the 
fundamental  through  the  higher-order  Kerr  response. 

The  above  qualitative  difference  between  the  predic¬ 
tions  of  the  old  and  the  new  models  makes  the  compar¬ 
ison  of  yields  of  the  third  and  fifth  harmonics  a  reliable 
test  of  the  new  theory  that  is  largely  insensitive  to  the 
variations  of  the  material  parameters  and  to  a  particular 
experimental  geometry. 

We  point  out  that  an  experiment  on  the  generation 
of  third  and  fifth  harmonics  in  an  Argon  cell  using  in¬ 
tense  pump  pulses  at  800  nm  wavelength  has  been  re¬ 
ported  much  before  the  HOKE  controversy  had  even 
emerged  [17].  The  results  of  that  work  have  been  re¬ 
cently  used  to  argue  in  favor  of  the  validity  of  HOKE  [18] . 
In  our  view,  [17]  is  not  straightforwardly  applicable  to 
the  resolution  of  the  HOKE  controversy  for  the  follow¬ 
ing  reasons.  Experiments  reported  in  [17]  were  certainly 
not  designed  to  support  or  disprove  the  HOKE  theory, 
which  was  nonexistent  at  the  time  when  [17]  was  pub¬ 
lished.  Instead,  that  work  targeted  the  demonstration  of 
highest  possible  yield  of  fifth  harmonic.  The  wavelength 
of  the  generated  fifth  harmonic  was  about  162  nm,  which 
is  in  vacuum  UV,  thus  making  its  reliable  and  quantita¬ 
tive  characterization  not  straightforward.  Furthermore, 
the  extended  interaction  length  and  the  associated  phase 
mismatch  between  the  pump  wavelength  and  the  gener¬ 
ated  harmonics,  together  with  the  fact  that  the  actu¬ 
ally  detected  portions  of  the  harmonics  were  not  exactly 
known,  made  the  interpretation  of  those  experiments  in 
the  context  of  the  resolution  of  the  HOKE  controversy 
not  reliable. 

The  experiment  reported  here  is  free  from  the  above 
complications.  In  our  case,  third  and  fifth  harmonics  are 
generated  using  millijoule- level  femtosecond  laser  pulses 
at  2.2  fi m  wavelength.  This  choice  of  pump  wavelength 
is  dictated  by  three  considerations.  First,  both  third  and 
fifth  harmonics  of  this  wavelength  fall  into  the  near  UV 
-  visible  spectral  range,  making  them  easily  detectable. 
Second,  the  phase-matching  distance  between  the  pump 
beam  and  both  of  the  harmonics  in  air  at  normal  at¬ 
mospheric  pressure,  at  this  pump  wavelength,  is  esti¬ 
mated  at  over  10  centimeters.  Thus  using  a  short  in¬ 
teraction  length  of  less  than  1  centimeter  in  our  ex¬ 
periment  eliminates  any  potential  issues  associated  with 
phase  mismatch  between  the  pump  and  the  harmonics. 
Third,  the  effect  of  the  higher-order  Kerr  terms,  accord¬ 
ing  to  [14],  should  become  even  more  pronounced  in  the 
mid-infrared,  compared  to  the  800  nm  wavelength  range 
commonly  used  in  the  high- intensity  femtosecond  laser 
experiments. 

Our  experimental  setup  is  schematically  shown  in 
Figure  1.  Femtosecond  pulses  at  2.2  /im  center  wave¬ 
length  are  produced  by  a  commercial  high-energy  Op¬ 
tical  Parametric  Amplifier  (OPA,  Palitra-HE  model  by 
Quantronix  Corporation).  The  OPA  is  pumped  by  about 
40  fs-long,  20  mJ  pulses  at  800  nm  wavelength  generated 
by  an  ultrafast  Ti: Sapphire  amplifier  chain  operating  at 
a  10  Hz  pulse  repetition  frequency.  The  OPA  is  tuned  to 


Fig.  1.  (Color  online)  Schematic  of  the  experimental 
setup. 


produce  signal  and  idler  pulses  at  1.26  /im  and  2.20  /im 
center  wavelengths,  respectively.  The  output  beam  from 
the  OPA  has  a  diameter  of  about  10  mm.  Both  signal  and 
idler  pulses  have  energies  of  about  1.5  mJ.  The  duration 
of  the  pulses  is  estimated  by  the  OPA  manufacturer  at 
between  40  and  60  fs  [19]. 

The  co-propagating  signal  and  idler  beams  are  focused 
by  a  lens  with  20  cm  focal  length.  Immediately  after  the 
lens,  the  1.26  fi m  signal  is  reflected  by  a  dichroic  beam¬ 
splitter  operated  at  near  normal  incidence.  This  dichroic 
(CVI  -  Melles  Griot  Corporation,  part  number  BBDS), 
is  a  dielectric  low-pass  filter.  At  normal  incidence,  it  has 
a  transmission  cutoff  at  the  wavelength  of  about  1.3  /im. 
The  idler  pulse  at  2.2  fi m  is  transmitted  by  this  mirror 
with  about  90%  efficiency,  while  the  transmission  for  the 
signal  beam  at  1.26  /im  is  less  than  0.1%.  For  all  wave¬ 
lengths  below  1  /im,  transmission  is  less  than  0.05%. 

The  energy  of  the  focused  2.2  /im  idler  pulse,  after  the 
dichroic,  is  measured  at  about  1.25  mJ.  It  is  important 
that  the  dichroic  is  oriented  with  its  substrate  facing 
the  incident  beam.  In  this  case,  the  visible  and  UV  light 
possibly  generated  inside  the  substrate  is  reflected  by  the 
coating,  and  only  the  2.2  /im  idler  pulse  is  passed  into  the 
interaction  zone. 

The  third  and  fifth  harmonics  of  the  2.2  /im  pump  are 
generated  in  the  ambient  air,  near  the  focal  plane  of  the 
focusing  lens.  The  estimated  intensity  of  the  pump  beam 
near  the  focus  is  above  the  turnover  point  for  the  alleged 
sign  reversal  of  the  Kerr  effect  in  the  HOKE  model  and 
also  above  the  ionization  threshold  of  air.  The  length 
of  the  short  and  faint  plasma  spark  generated  in  the 
vicinity  of  the  interaction  zone  is  less  than  1  cm.  Under 
these  conditions,  the  ratio  of  the  third  and  fifth  harmonic 
yields,  according  to  the  HOKE  theory,  should  approach 
one. 

The  third  and  fifth  harmonics  of  the  pump  are  re¬ 
flected  by  another  dichroic  mirror,  which  is  identical  to 
the  one  used  to  block  the  signal  pulse  at  1.26  /im  from 
the  OPA.  The  second  dichroic  is  operated  with  the  coat¬ 
ing  facing  the  incident  beam  and  at  a  45°  angle  of  inci¬ 
dence.  This  mirror  efficiently  reflects  the  generated  har¬ 
monics,  but  blocks  over  90%  of  the  2.2  fi m  pump.  The 
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Fig.  2.  Normalized  spectra  of  third  harmonic  (left)  and 
fifth  harmonic  (right)  generated  by  tightly  focused  fem¬ 
tosecond  pulses  at  2.2  /im  wavelength  in  air,  on  linear 
scale. 


residual  pump  is  further  attenuated  by  about  five  orders 
of  magnitude  by  a  3  mm -thick  color-glass  filter  plate 
(KG5  filter  by  Schott  Glass  Corporation)  placed  after 
the  dichroic. 

After  collimation  with  another  lens,  the  signals  at 
third  and  fifth  harmonics  are  spatially  separated  by  a 
300  lines  per  millimeter  transmissive  diffraction  grating. 
The  spectra  of  the  two  harmonics  are  recorded  by  a  fiber- 
coupled  spectrometer  (model  USB4000  by  Ocean  Optics 
Corporation),  with  the  aid  of  two  additional  lenses  that 
separately  focus  the  harmonic  beams  after  they  are  sep¬ 
arated  by  the  grating.  These  spectra  are  shown  in  Figure 

2.  They  are  centered  at  733  nm  and  438  nm,  very  close 
to  one-third  and  one-fifth  of  the  pump  wavelength  of 
2.20  /im. 

The  emission  bandwidths  of  the  detected  third  and 
fifth  harmonics  are  about  20  nm  and  11  nm,  respectively. 
Assuming  that  the  harmonic  pulses  are  both  Gaussian 
and  close  to  transform  limit,  the  durations  of  the  two 
pulses  can  be  estimated  as  40  fs  and  25  fs.  The  ratio  of 
these  pulse  durations  is  in  the  ballpark  of  yj 5/3,  which 
would  be  expected  in  the  case  of  the  perfectly  phase- 
matched  harmonic  generation  with  non-depleted  pump. 

The  energies  of  the  two  harmonic  pulses  are  meas¬ 
ured  by  dedicated  silicon  detectors  that  were  calibrated 
against  a  pyroelectric  detector,  using  reference  laser 
pulses  at  800  nm  wavelength.  Accounting  for  all  losses 
experienced  by  the  generated  harmonics  in  our  setup, 
the  energies  of  the  third  and  fifth  harmonics  are  found 
equal  to  1.9 /iJ  and  0.4  nJ,  respectively.  These  pulse  en¬ 
ergies  correspond  to  yields  of  about  1.5  •  10-3  for  third 
harmonic  and  3.0  •  10-7  for  fifth  harmonic,  with  respect 
to  the  energy  of  the  pump  pulse  at  2.2  /i m  wavelength 
in  the  interaction  zone.  The  ratio  of  yields  for  the  two 
harmonics  is  about  2  •  10-4,  in  contradiction  with  the 
HOKE  theory. 

To  verify  that  the  generated  harmonics  were  in  fact 
originating  from  the  tightly  focused  interaction  zone  and 
not  from  any  optics  in  the  setup,  we  substituted  the  fo¬ 
cusing  lens  with  a  glass  plate  with  the  same  thickness  as 
the  lens.  In  that  case,  the  intensities  of  both  third  and 


fifth  harmonics  of  the  pump  fell  below  the  detectability 
limit. 

In  conclusion,  we  have  conducted  a  dedicated  experi¬ 
ment  in  order  to  validate  or  disprove  the  existence  of  the 
higher-order  Kerr  effect  proposed  in  [1].  Our  test  was 
based  on  the  qualitative  difference  that  the  inclusion  of 
the  higher-order  Kerr  terms  has  on  the  yields  of  the  third 
and  fifth  harmonics  generated  by  tightly  focused  ultrain¬ 
tense  laser  pulses  at  2.2 /rm  wavelength  in  ambient  air. 
Our  results  unambiguously  disprove  the  existence  of  the 
high-order  Kerr  response,  even  in  the  transient  regime. 
Our  conclusion  is  in  agreement  with  another  recently  re¬ 
ported  dedicated  test  of  HOKE  [15]. 
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